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Abstract 

Fit;o) is a Mossbauer fitting and analysis program written in Borland Delphi. It 
has a complete graphical user interface that allows all actions to be carried out via 
mouse clicks or key shortcut operations in a WYSIWYG fashion. The program does 
not perform complete transmission integrals, and will therefore not be suited for a 
complete analysis of all types of Mossbauer spectra and e.g. low temperature spec- 
tra of ferrous silicates. Instead, the program is intended for application on complex 
spectra resulting from typical mineral samples, in which many phases and differ- 
ent crystallite sizes are often present at the same time. The program provides the 
opportunity to fit the spectra with Gaussian, Lorentzian, Split-Lorentzian, Pseudo- 
Voigt, Pseudo-Lorentz and Pearson- VII line profiles for individual components of 
the spectra. This feature is particularly useful when the sample contains compo- 
nents, that are affected by effects of either relaxation or interaction among particles. 
Fitted spectra may be printed, fits saved, data files exported for graph creation in 
other programs, and analysis tables and reports may be exported as plain text or 
I^TgX files. With Fit;o) even an inexperienced user will soon be able to analyze and 
fit relatively complex Mossbauer spectra of mineralogical samples quickly without 
programming knowledge. 
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1 Introduction 



During work with analysis of Mossbauer spectra we found that existing anal- 
ysis software packages however competent (i.e. Mfit [1], Recoil [2], MacFit [3], 
MossWinn[4]) could not fulfill our needs of being simple and yet very flex- 
ible. For fitting of spectra with many components, as is typical for samples 
with complex mineral assemblies, advanced programs are in many cases too 
detailed to be really useful. 

We therefore developed this program package which attempts to mitigate these 
issues, and thus provides a powerful tool for rapid fitting of complex Mossbauer 
spectra. 

The program has been downloaded around 100 times in several versions, and 
an estimate of regular users is around 20. The program package has been used 
in the preparation of several articles or talks [5,6,7,8,9]. 



2 Program summary 

Fit;o) is a program for fitting and analyzing transmission and scattering ge- 
ometry ^^Fe Mossbauer spectra of metals, alloys and mixtures of ferric oxides, 
oxyhydrates and silicates. The program can be expanded to handle other Moss- 
bauer isotopes as well. It has a complete graphical user interface, which allows 
all actions to be performed via mouse clicks or key shortcuts. The program 
accepts currently both particular plain text files (.exp-files, see [10,11] for a 
definition) and comma-separated files e.g. counts vs. velocity (velocity;counts) 
as data input. It is the intention to release the source code under an open li- 
cense. 

Mossbauer spectra can be fitted using singlets, doublets and sextets (all with 
a selection of line profiles of which Lorentzian is default) and the fit model 
can be saved for later reloading for example for further refinement of the fit. 
Fitted and unfitted spectra can be saved as reports and saved as both plain 
text and in customizable I^T^]Xformat, and spectra can be printed with a fit 
report. 

The program contains a precompiled list of Mossbauer data of many common 
iron compounds with an option for the user to edit or add new compounds. 

Calibration data for Mossbauer experimental setups can be extracted from 
calibration spectra, and saved in separate calibration .cal-files (see [10,11] for 
a description). 
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The settings window contains a wide range of customizable settings through 
which the user may customize the appearance and behavior of the program. 
Backwards compatibility is ensured, as file formats are kept unchanged and 
user-changed files are not overwritten. 

The web update option provides an easy and automated way to keep the 
software updated with the latest program version. 

The program does not perform complete transmission integrals, and will there- 
fore not be suited for a complete analysis of all types of Mossbauer spectra. 
Also the fitting does not use the complete spin-Hamiltonian (but approxima- 
tions) so for example magnetic (low-temperature) spectra of silicates cannot 
be analyzed by the programme. 

Some of Fit;o)'s main features are: 

• Microsoft Windows 2000/XP compatible. 

• Easy to install, no external dependencies, and safe uninstallation/removal. 

• Complete point-and-click graphical user interface. 

• Easy saving and loading of fit models. 

• Most program parameters as customizable. 

• Working with several spectra (mdQ at a time is possible. 

• Created with object-oriented programming. 



3 Description 

3. 1 Prerequisites 

As noted in the introduction, execution of Fit;o) requires an updated Microsoft 
Windows 2000, Microsoft Windows XP operating system. A screen resolution 
of a least 1024x768 with a 16 bit color depth is recommended. 

3.2 Installation and execution 

The program is available at http://hjollum.com/jari/zzbug/fit/. After 
installation and execution the main options are: 

• Open the selected type spectrum by selecting an appropriate icon in File 
menu. 

^ Multi Document Interface 
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• To open a file in text mode, use the Open text file option and to view 
the application log, use the Open log option. 

A more detailed description of the program interface can be found in the 
manual [10] or the web site [11]. 



3.3 The layout of the graphical user interface 



After the first execution of the program the user is presented to the graphical 
user interface shown in figure 1. The main window of Fit;o) is divided into five 
areas (listed from top down), the menu bar, the tool bar, the work area, the 
window panel and the status panel. 



34 Fitting 




[■^ Fitting :A1 DDD2.EXP 



|open Mosshauer transmis5iein spectrum ^ 

Fig. 1. Opening a spectrum file. 

When choosing to open a spectrum for fitting, the open spectrum form (figure 
2) appears, in which the spectrum to be loaded is chosen. After choosing a 
spectrum the fit form (figure 3) is opened and the spectrum is loaded. 

As fitting components, it is possible to insert either simple (fundamental) 
components (singlets, doublets, sextets), or to insert predefined mineral com- 
ponents. 

Simple components are inserted by clicking one of the spectrum icons in the 
toolbar, and then marking, position, intensity and width of the components 
on the graph. The parameter that the user is expected to mark is indicated 
at the cursor. 
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Fig. 2. Choosing the spectrum to open. 




■t^|Filliiig'Aia002 EXP 



Fig. 3. The fit form with a loaded spectrum. 

Predefined (but editable) mineral components are inserted by clicking the 
cylinder icon in the toolbar, and choosing the mineral. After choosing, one 
has to mark the intensity of the component on the graph using the mouse 
cursor. 

The fitting algorithm is chosen in the Fit Controls combobox in the lower 
left corner of the form. The fitting process is started and aborted by using the 
adjacent start and stop. 



3.5 Calibration 



Fitting components are inserted by clicking the button in the toolbar with the 
sextet icon (see figure 5). When the components have been inserted, either a 
single fitting run or a series of repetitive fitting runs can be started, through 
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Fig. 4. The fit form with an inserted component. 
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Fig. 5. T/ie calibration form before inserting the calibration fit components. The insert 
button is marked. 





Fig. 6. The calibration form after insertion of the calibration fit components. 

the start icons. The result can be saved and exported by chcking the 'notepad' 
icon. 
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Fig. 7. The minerals administration window. 

3. 6 Minerals 



The program contains a list of common and well-known minerals (figure 7) , 
which can be used for as starting point fitting. The mineral data are from [12]. 
It is possible to edit the listed minerals and save these modifications, and to 
enter new minerals. This can be done either through the New, Edit or Delete 
options, or by editing the minerals file manually. The name and location of 
the minerals file is displayed in top of the window. 



4 Project planning 



In the following we discuss the properties of the program as used for analysis 
of a transmission/absorption spectrum. However, the discussion is completely 
valid also for scattering/emission spectroscopy, since there will be only few 
details differentiating these. 



4-1 Line shapes 



Ideally each Mossbauer absorption line has a Lorentzian line shape, (see table 
1). However, there are physical effects that might disrupt or distort the ideal 
case. For these cases other line shapes are necessary. 

In practice most absorbtion lines are Lorentzians, which may be sligthly 
smeared by a Gaussian due to for instance temperature fiuctuations, leading 
to the Voigt line shape. The Voigt line shape is a convolution of a Lorentzian 
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and a Gaussian and can be expressed as [13] 



'Voigt 



■'Gau 



'v')LLor{v - V')dv'. 



(1) 



However, this can not be solved analytically, and therefore several approxima- 
tions to this profile exist. We have implemented three of these: Pseudo-Voigt, 
Pseudo-Lorentz and Pearson- VII, which are listed in table 1. Especially the 
Pearson- VII line shape is difficult to implement, since it involves the imple- 
mentation of the Gamma function [14], which was implemented using Stirlings 
approximation [15] 

TsUrUng{z) ^ ^ siuh ^ + ^LT) , (2) 

which is accurate to a least 5 digits for z>l. 

j,^(js Split-Lorentz line profile 
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X 

Fig. 8. An example of a spectrum with Split-Lorentzian line profiles. The outer slope 
of an absorption line is defined as the slope farthest from the center of the component. 
This of course only applies to doublets and sextets. 

The Split-Lorentzian (see figure 8) line shape is commonly used for samples 
in which small-particle effects or interactions distort the line shapes in such a 
way that the outer slope is steeper than the inner slope. The Split-Lorentzian 
line shape is made up of the outer side using one Lorentzian line shape, and 
the inner side using another shape. The amplitude and center of the two line 
shapes are kept equal, but the widths vary. The relationship between the 
widths is usually controlled by defining a width parameter and an asymmetry 
parameter. The asymmetry parameter is usually known as b. 

When using the Split-Lorentzian line shape in doublets and sextet (see section 
4.2) the line shapes are mirrored with respect to the center of the component, 
and the left side of a line shape is identical to the right side of its mirrored 
twin. 




8 



The line shapes implemented for the Mossbauer analysis and the properties 
of these are presented in table 1. 
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Table 1 

Line shapes used for Mossbauer analysis. Reproduced from [16]. 



4-2 Fit components 

The three fundamental Mossbauer fitting components are the singlet, the dou- 
blet and the sextet with 1, 2 and 6 absorption lines, respectively. All of these 
components are characterized by their isomer shift (IS). 

The splitting of the doublet depend on the quadrupole interaction, and is in 
the doublet quadrupole splitting (QS). 

The singlet, doublet and the sextet have furthermore the line shape a parame- 
ter (6, 77, a and m in table 1) in common. The Split-Lorentzian, Pseudo-Voigt, 
Pseudo-Lorentz and Pearson-VII line shapes use this property. 

The magnetic Zeeman interaction, the magnetic hyperfine field (HF), is only 
used in the sextet, and needs only to be implemented in the sextet. For sextets 
this software - presently - can handle only spectra in which the magnetic 
Zeeman interaction is dominant, so that the quadrupole interaction can be 
treated as a small perturbation on the magnetic interaction. In the case of the 
sextet QS indicates the quadrupole shift in contrast to the doublet quadrupole 
split. 

However, since the properties HF, and QS, are the only properties which do 
not apply to all components, we have provided an interface to them in all 
components, but disabled them where they are not needed. As mentioned, we 
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do not use the complete spin-Hamiltonian for calculations of the line positions 
in the, but rather approximate expressions, which however, have been quite 
sufficient in most all cases. 

The singlet consists of a single line, and the position of the single line is the 
same as the isomer shift. The expression for the singlet is given by 

yi = L{vi-vo), (3) 



where L is the line shape function providing the value of Hi as a function of 
Vi, and Vq is the value of the isomer shift. 

The doublet consists of two lines, and the centers of the lines are placed at a 
distance of QS apart, around the isomer shift. The expression for the doublet 
is given by 

Vi^ L{vi-VQ-^) + L{vi-vo + ^), (4) 
where vqs is the value of the quadrupole splitting. 

The sextet consists of six lines, and the centers of the lines are placed at 
positions dictated by the isomer shift, quadrupole shift and hyperfine field. 

The expression for the sextet is given by 

y, = Y.Hvi -Vo + + i-lf+'-^^^k,,^ ■ Bhf), (5) 

i=i ^ 



where bmnki is a delta function, which is 1 when i has one of its values, kq^w 
is a Mossbauer proportionality factor, composed of several physical constants. 
The three values that apply are k\(^, /c2,5, /c3,4 for the paired lines. B^p is the 
hyperfine field. The /c, j values are calculated from the general expression for 
the sextet energy levels, which is given by 

Em = -gPnmB, (6) 



where g is the Lande factor, /3„ is the nuclear magneton, and m is the m- 
quantum number. The energy splitting is calculated by 

|AE(ij)| = \AEi\ + \AEj\ = 2\AEi\, (7) 

where = (1, 6), (2, 5), (3, 4). Applied to the line pairs (1,6), (2,5), (3,4) 
then energy splitting is calculated from 
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|A%,)| = -2(/3n5(^e^ - 9g\)) = PnB ■ Qij, (8) 

where Qg = 0.181208 and ge = —0.10355 are the g-factors for the ground state 
and the excited state respectively, values from [12], and the gij factors. 



516 = -2(5e^ - gg^) = 0.491858, (9) 
525 = -2(ye^-54) = 0-284758, 

934 = -2(5e^ - 9gl) = 0.077658, (11) 

(12) 

are g-factors extracted from the g-factor for the excited and ground state 
and £'0 is the transition energy from 7 = | to / = | for ^^Fe Mossbauer 
spectroscopy. 

The actual calculation being performed is 

Q 

VhF = gijPnBHF -TTET: (13) 



using unit conversion this becomes 

VHF = 9ijBHF^^^^^ (14) 

5.0505 • 10-2^ J/T • 3.0 • 10^1 mm/s 
- Qij^HF 2 . 41 j^eY . ^ 5322 • lO-i^J/eV ^ ^ 

= QijBHF ■ 0.32794 mm/sT, (16) 
where e is the elementary charge. Now the fc^^^'s can be calculated as 

fci,6 = 1.61299- 10"^ mm/sT, (17) 
A;2 5 = 9.33835- 10"^ mm/sT, (18) 
A;3,4 = 2.54672 - 10"^ mm/sT. (19) 

4.3 Fitting spectra 



The program will analyze Mossbauer spectra by fitting a model set by the 
user, and it will report the result of the fit. 
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Before any analysis can commence, a background level has to be established/calculated. 
The background level is calculated as the mean of the value of the outermost 
8 channels on each side of the spectrum. Transmission spectra contain absorp- 
tion lines, which have fewer counts than the background level. The expression 
for components provided above has to be subtracted from the background to 
produce a model data series. 

The program can also be used for analysis of reflection spectra. Reflection spec- 
tra also contain a background level, to which the reflection lines are added. To 
produce a fltting model, the fltting components are added to the background 
level. 

4-4 Calibration spectra 

The program derives the calibration parameters of the spectrometer from au- 
tomatic analysis of calibration spectra. The found calibration parameters are 
then used in the analysis of spectra obtained under circumstances identical 
to those of the calibration spectrum. In this program version (1.0.0.63) only 
calibration of linear velocity profiles are implemented. We are well aware that 
sinusoidal velocity profiles are common in the Mossbauer community, and a 
later version may be adapted for spectra obtained in the mode. 

Contrary to a normal spectrum for Mossbauer analysis, a calibration spectrum 
is analyzed unfolded with the channel numbers used as reference. 

The calibration spectrum usually consists of the spectrum of a thin iron (a-Fe) 
foil (12.5 /im) at 295 K. But other calibration materials can also be used, such 
as stainless steel and other well characterized iron compounds. 

The unfolded calibration spectrum consists of a background level, and 12 
absorption lines on the background level figure 5. The calibration is used 
for finding 3 unknown parameters of the experimental setup, and to provide 
information on any imperfections of the source/drive system i.e. increased line 
widths, nonlinearity etc. The three unknown parameters are listed below: 

The folding channel is the channel (integer) around which the data chan- 
nels are folded. Operationally (for historic reasons and backwards compati- 
bility) it is defined as the channel number, which is added to channel number 
1, e.g. it is the double of the actual folding channel. In a 512 channel setup, 
it has typically a value of 510-514. It is found by calculating the mean of 
positions of the lines (1,12), (2,11), etc. 

The zero velocity channel is defined as the detector channel position (dec- 
imal, 5 digits), of the symmetry center of the 295 K spectrum of a thin foil 
of CK-Fe. It is found in each of the spectrum halves separately. It is found 
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by locating the mean center of the lines (1,6), (2,5) and (3,4) pairwise, and 
correspondingly for the second half of the spectrum, and then finding the 
common center. 

The calibration constant is the link between the channel number and the 
source velocity. It is found as the mean value of the known hyperfine field 
of the absorber material {Bhf = 33.02 T for a-Fe at 295 K) divided by the 
distances between the transition pairs, = (1, 6), (2, 5), (3, 4), and scaled 
by the Mossbauer constants (17). The exact formula is 

6 Eq ■ C \Pj — Pi Pj+6 — Pi+6 / 

where Pi is the channel position of the i'th absorption line. 

The background line has a shape that is dependent on the mode of operation 
of the spectrometer (e.g. frequency and velocity range) and geometry of the 
experimental setup. The radioactive source is placed inside a collimator, and 
since the source is oscillating with near constant acceleration, the solid angle 
seen by the source will vary as a function of position. The radiation detected 
is proportional to the angular area. The baseline shape is approx. that of two 
parabolas in succession with positive and negative a values: 

A = +a ■ {v - - b{v - Cm) + c, I <v < cm, (21) 
A = —a ■ {v — cuf' — b{v — Cm) + c, cm < v < N, (22) 

where A is the angular area, cm is the mirror channel, v is the channel number, 
N is the number of channels, and a, b and c are the parabola parameters. 

Besides the above calculations, the calibration is able to find the 12 absorption 
lines (''"^Fe), and fit these and the baseline to find the calibration constants. 

As the equipments available during the development only used a linear velocity 
drives, only this case has been implemented in the program, but this is easily 
expanded, when the apparatus data are available. 



5 The program structure 

The program is built entirely on object oriented technology (OOT). There are 
several advantages in using OOT: 

• Inheritance makes it easy to create a child class, which inherits most of 
its properties and methods from its ancestor, but introduces some new, or 
changes the internal behavior in some areas. 
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• Isolation makes it easy to correct errors or undesirable behavior, without 
affecting other parts of the program, thereby minimizing program errors. 
Furthermore it makes it easy to extend functionality. 

There are also disadvantages in choosing OOT: 

• Processing speed will in most cases be lower than that of a procedural pro- 
gram since the overhead will be larger, and the amounts of data manipulated 
in most cases will be larger. 

• The implementation process will take longer and the source code will be 
larger, since similar behavior may have to be implemented several times in 
order to avoid code shared between classes. 

The program was implemented using Borland Delphi 5, 6 and 7, and should 
be executed on a PC running Microsoft Windows 2000 or Microsoft Windows 
XP. During implementation the built-in classes of Delphi have been used as 
programming base. For increased flexibility the program is designed to be a 
multiple document interface (MDI) program, meaning that it would be possi- 
ble to work with several 'documents' (spectra) simultaneously. 

5. 1 The fitting form 

The Mossbauer fitting form is implemented as the class TdfmFitForm (see 
figure 9). It is the graphic user interface, through which the user makes input to 
and receives output from the spectrum analysis. It contains two very important 
classes TGraph and TFitGraph. 

After the FitForm has been created and initialized, it loads a data-file (which 
the user has chosen), and calculates the Mossbauer spectrum background level. 
The spectrum is folded during loading. 

Whenever a change has been made to one of the fitting components, or a new 
fitting component has been inserted, several calculations have to be performed: 

• The sum data, SumData, which contains the sum of the fitting components 
is calculated as 

Ntet's 

Si^T. Tj,,, (23) 

3=1 

where Si is the sum at the i'th point, Tj i is the value of the j'th fit component 
(singlet, doublet, sextet), at the i'th point. 

• The rest data, RestData, which contains the difference between the spec- 
trum data and the sum is calculated as 
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Fig. 9. The logical/workflow diagram of the TdfmFitForm class, used for fitting Moss- 
hauer spectra. The upper part is the logical/workflow diagram, and the symbol expla- 
nations are shown in the lower part. 



Ri = Di — {Bi — Si) 
R^ = D,-{B, + S,) 



Di — Mi for transmission spectra, 
Di — Mi for scattering spectra, 



(24) 
(25) 



where Ri is the rest at the i'th point, Di is the i'th data point of the spectrum 
data, Bi is the i'th point in the baseline and Mi is the model data at the 
i'th point. 

The ChiSqr, and which are indicators for difference between the 
model and the data, and used for fitting, is calculated. The misfit is also 
given as an output parameter. It can be summarized mathematically to 



X 



^valid p2 



'^valid i=\ -^i 



(26) 



since Ei = ^/Di and riband is the number of valid channels in the spectrum. 
Furthermore we calculate the reduced chi-squared by 



1 'Invalid p2 

/ y 7-) ) 



(27) 
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where nfree is the number of free fit parameters. As a goodness of fit pa- 
rameter we use the misfit m [17] defined by 



m 



n-valid (Bj—Dj) 



i = l 



(28) 



Two fitting algorithms, random walk [16] and amoebe [18], have been imple- 
mented in the fitting form. First we will go through the random walk algo- 
rithm. 



ChiSqrBest:=ChiSqr 



Save initial settings. 



Outer loop variable 



lnr)er loop variable 



Generate random fit parameters 



Calculate values 



If (ChiSqr/ChiSqrBest)<0.999 

\ 

t 

ChiSqrBest;=ChiSqr 

[^ KeepRandomFit ^ 

1 H>--inner 



User termination 



Evaluate parameters 



Saves the parameters of 
the best fit 



Fig. 10. A graphical illustration of the RandomFit fitting algorithm. The steps marked 
with octagonal boxes indicate that this step iterates through all TFitTet's in the fit. 

The random walk algorithm is a variant of the "simulated annealing"-algorithm 
used in many areas of physics. It differs however, in that the tunneling possi- 
bility has been eliminated. The RandomFit algorithm is presented graphically 
in figure 10. The algorithm is based on two loops, the outer and the inner. 
The core of the algorithm is thus run 



N 



n. 



outer 



Til 



(29) 



times. The algorithm starts with setting the variable step to 1. Each time the 
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outer loop is traversed once, step is set to one half of its previous value. The 
philosophy is to find the hitherto best fit for each time the outer loop has been 
traversed, and thereafter to narrow the interval in which the parameters are 
allowed to vary, before traversing to the outer loop again. 

The algorithm is based on the assumption that the best fit after n,„„e,. guesses, 
lies in the vicinity of the best fit. The fitting parameters all have an interval 
to which they are restricted during fitting. Using the result of the sum 

oo 

^2-^ = 2, (30) 

1=0 



this can be maintained/fulfilled. If A is the fitting variation interval of a fitting 
parameter, the following formula 

£ A-2-*<2-A (31) 

1=0 



provides a way to implement our version of the random walk algorithm without 
violating the restrictions: 

J2 1/2- A-2-^ < 1/2-2- A = A. (32) 

i=0 

Before the first traversation of the inner loop, the variation intervals are there- 
fore set to half of the original, corresponding to the formula above. This guar- 
antees that the fit parameters are kept within the allowed fit interval. 

The amoebe algorithm [18] works in many areas much like the random-walk 
algorithm. It has an outer loop, which repeats riouter times, or until the im- 
provements are no longer significant. The inner loop basically works same way 
as in random walk. 

Besides the evaluation of the improvements mentioned, the difference between 
the random walk and amoebe is that after the best fit has been found using 
the inner loop, some preprogrammed variations of the best fit are tried. For 
all paramters the following is tried individually: Try to the fit parameter value 
in the opposite direction, with the same amount. If this gives a better fit, 
it tries to shift the parameter value a step further in the same direction, to 
see if it should give an even better fit. This tests if the point of origin was a 
local maximum, and if the fit has found the worst local minimum. If so, the 
algorithms tries if a further step in the same direction gives a better fit. 

If the opposite step did not yield a better result, the algorithm tests to see if 
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Fig. 11. A graphical illustration of the preprogrammed steps in the amoebe fitting 
algorithm. 

it has stepped too far to find the currently best fit, and tries to see if there 
is a better point halfway between the origin and the currently best fit. If all 
of these attempt do not result in a better fit the step variable is divided by 
2, and the content of the outer loop is traversed again. If a better fit is found 
step is not changed. The preprogrammed steps are illustrated in figure 11. 

The outer loop will terminate when it has been traversed 20 times, terminated 
by the user or when the rtol variable is larger than the ftol. ftol is a 
measure of the minimum relative improvement required, for improvements to 
be considered significant, rtol is a measure of the improvement of the best 
fit compared to the worst fit in the current traversation of outer loop, rtol is 
calculated from: 



rtol 



(33) 



The source code is presented graphically in figure 12. 



5.2 The calibration form 



The Mossbauer calibration form is the graphic user interface, through which 
the user makes input to and receives output from the calibration process. 
Generally it resembles the fitting form in many points, and the description 
of common features will not be repeated. The data files are, however, loaded 
unfolded. 

Whenever a change has been made to one of the fitting parameters, the same 
calculations as in the fitting form have to be performed. SumData, PlotSumData, 
RestData and PlotRestData, are calculated in the same way. The only diff'er- 
ence with regards to the fitting form is the variable baseline in the calibrations 
form, which has to be taken into account. 
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start 

i_ 

Initialize 



iter:=0 



31 



ChiSqrNWorst:=ChiSqrBest 
ChiSqrWorst:=ChiSqrBest 

i 



CalGVaiues(step) 

i 



if ChiSqr<CliiSqrBest 



if CfilSqr>ClniSqrWorst 



if ChiSqr>CliiSqrN Worst 



if ChiSqr<CliiSqrBest 



if ChiSqr>CliiSqrNWorst 



if ChiSqr<ChiSqrBest 






1 

Tme 
▼ 


KeepBestFit 







if ChiSqr<ClniSqrBest 



iter:=iter+1 
if (notdone)and(iter<=itmax) 



End 



Fig. 12. A graphical illustration of the AmoebeFit fitting algorithm. The steps marked 
with octagonal boxes indicate that this step iterates through all TFitTet's in the fit. 

Before the calibration process can commence, the fitting components have to 
be inserted, and the baseline parameters defined. The baseline is initially fixed, 
but is allowed to vary during the calibration. The 12 components are inserted 
at the 12 local extrema with the largest 'magnitude'. 

The local minima are found by traversing through the data points of the data, 
and for each data point it checks whether the current point is a local extremum 
by comparing it to its neighboring points. If the point is a local extremum the 
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TcHmCuslomCalSingletForm 



TdftnCalSingletForm 



TdfmCalFjlForm 



TFitCalBaseLine 



TCalBaseLme 



RawData j*- 
-t- j Sum Data ~| 
-■j FlotSumPala |^- 




A and B have references to, and use each olher. 



A has a reference to, and u; 



n ancestor of B and C. and references to A 
e to implemented as instances of B or C. 



B and C are instances of TType, with different roles. 



TTypes is a cciiection/iist/array of objects, if it references a class class, it contains arbifranly 
many objects of the class. If it references an object, the object is a fixed member of the 
collection. 



Fig. 13. The logical/workflow diagram of the TdfmCalForm class, used for fitting 
Mossbauer spectra. The upper part is the logical/workflow diagram, and the symbol 
explanations are shown in the lower part. 



magnitude is calculated as 



yi+2 + Vi+i + Vi-i + yi-2 



pivot 



(34) 



where i/i is the i'th point. After finding the all local extrema, the list of extrema 
is sorted, so the local extrema with the largest magnitude are placed first in 
the list, and the predefined number of largest local extrema are returned, and 
used for inserting the fitting components. 



The random walk in the calibration form is very similar to that of the fitting 
form. The only difference lies in that the inner loop fits the components pa- 
rameters when the inner loop counter j is even and the baseline parameters 
when it is odd. This is done since experiments during programming shows that 
this improves the intelligence and speed of the calibration procedure. 
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6 Conclusion 

We have presented the functionality and the basic principles behind the Fit;o) 
program. We believe that by the program we have made available a valuable 
tool for general simple Mossbauer analysis of complex samples. 
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8 Additional information 

The program is available from the Fit;o) homepage at [11]. There is a forum 
section, a FAQ, as well as a short manual available. 
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A Examples 



There are a few .exp-files included in the installation package. These are test 
spectra, which can be used for practicing. A description of the contents of the 
files is given below: 

• AlOOOO.exp, AlOOOOel . exp, AlOOOOer.exp are simulated/generated test- 
ing spectra. They contain a perfect a-Fe sextet on both, right and left side 
respectively. 

• AlOOOl.exp is a spectrum of soil from Salten Skov, Arhus, Denmark. 

• A10002.exp is an a-Fe transmission calibration spectrum from a Fe foil. 

• BlOOOl.exp is an a-Fe scattering calibration spectrum from a Fe foil. 

• B10105A.exp BlOOOl.exp is an a-Fe scattering calibration spectrum from 
a 25 ;U Fe foil, with an perpendicularly applied magnetic field. 

• B10106a. exp is an a-Fe scattering calibration spectrum from a 25 /i Fe foil. 

• HA668.exp, HA679.exp, HA869.exp, HA900.exp, HBSOO.exp, HB506.exp, 
HB522.exp and HE0195.exp are real spectra of various samples, both natu- 
ral and synthesized. 
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